Multi-object Spectroscopy of Young Star Clusters in NGC 4676 
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ABSTRACT 

Galaxy interactions are known to trigger starbursts. The young star clusters formed in merg- 
ers may be young globular clusters. The ages of these young star clusters yield the timing of 
interaction-triggered star formation and provide an important way to reconstruct the history of 
merging galaxies. Here we present the first results from our investigation into age and metallicity 
of twelve young clusters in the merging galaxy pair NGC 4676, using spectra from the multi- 
object Low Resolution Imaging Spectrometer (LRIS) on Keck. For ten clusters, comparison of 
the Balmer emission lines with model equivalent widths (EWs) yields ages less than 10 Myr. Two 
spectra display Balmer absorption lines typical of star clusters dominated by A-type stars, with 
estimated ages of about 170 Myr. These ages are comparable to the dynamical age of the tidal 
tails and are consistent with star formation triggered during the first passage of the pair. The 
locations of these two clusters in the tidal tails are generally consistent with predictions of shock- 
induced star formation models. One of these older objects appears unresolved on the image and 
is luminous enough to qualify as a young globular cluster. Using EWs of the diagnostic lines 
[O II] and [O III], we obtain oxygen abundances in the range 7.3 < 12 + log(0/H) < 9.0. These 
values show a nearly flat distribution along the northern tail, suggesting efficient gas mixing in 
the tail. 



Subject headings: galaxies: star clusters 
galaxies: interactions 

INTRODUCTION 



galaxies: individual(NGC 4676) — galaxies: abundances 



Gas- rich merging galaxies, e.g. NGC 4 038/39 
dWhitmore et alJl999L etc ) and NGC 3256 (|Zepf et al 
Il999t iTrancho et al.ll2006l ) , are now known to host 
numerous young massive star clusters that formed 
during the merger. These young star clusters 
formed during an interaction-induced starburst 
can provide a powerful diagnostic for t racing star 
formation history (e.g. ISchweizerlll982n . With the 
Mice as a starting point, our goal is to use multi- 



object spectroscopy to map the ages of young 
clusters throughout the bodies and tails of inter- 
acting galaxies, and thereby trace the history of 
star formation during the encounter. 

The two galaxies making up NGC 4676 were 
dubbe d "the Playing Mice" bv lVorontsov-Vervaminov 
(|1958l ). The pair is located at a,i2ooo.o = 
12H3 m 44 s , (5j2ooo.o = +31°00T8" femes et al 
200lh and has a recession velocity cz Q 
Adopting a Hubble constant of H 



1 The data presented herein were obtained at the W.M. 
Keck Observatory, which is operated as a scientific part- 
nership among the California Institute of Technology, the 
University of California and the National Aeronautics and 
Space Administration. The Observatory was made possi- 
ble by the generous financial support of the W.M. Keck 
Foundation. 



eeiokms- 1 . 

72 kms^Mpc -1 
places the Mice at a distance of 88.8 Mpc, corre- 
sponding to a scale of 1" = 0.43 kpc. A dynamical 
simulation suggested that the first pericentric pas- 
sage of the pair occurred about 170 Myr ago, and 
that the extended star formation observed in the 
system is b etter describe d by shock-triggered star 
formation ( Barnes! 20041 ). To test the validity of 
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this model, we obtained spectra of twelve young 
star clusters and associations in the Mice. Here 
we present the spectra and first results on the ages 
and metallicities of these clusters. 

2. OBSERVATION AND REDUCTION 

Observa tions were performed on 2005 May 4, 
with LRIS (|Oke et al.lll995l ) on Keck I. The target 
star clus ters were chose n from the brighter candi- 
dates in Ide Griis et al. (2003); we attempted to 
sample clusters from every part of the color-color 
diagram (see their Figure 4). T he astrometric so- 
lution in de Griis et al.l ( 2003h was based on the 
Hubble Space Telescope (HST) images, producing 
an intrinsic scale error of about 1% in the tabu- 
lated coordinates of these clusters. Corrections to 
the coordi nates were made u sing the USNO-A2.0 
catalogue ( Monet et all 19981 ) as an absolute astro- 
metric standard. Table 1 shows the corrected co- 
ordinates and measured properties of each cluster. 
A g-band image of the Mice taken by the Adva nced 
Camera for Surveys (ACS; iFord et al.l Il998h and 
the position of slits and clusters in the Mice are 
shown in Figure 1. Each slit width is 1", with 
lengths ranging from 5" to 16". Two slits captured 
spectra of additional off-target stellar populations 
and star clusters. The blue spectra were obtained 
with a dispersion of 1.09 A per pixel, a spectral res- 
olution of 5.95 A, and an approximate wavelength 
coverage of 3200 A to 5600 A. Nine 30 minute ex- 
posures were taken, together with standard star 
observations and short calibration exposures. Five 
of the exposures were partly obstructed by passing 
cirrus, thus only four frames were used to obtain 
the spectra analysis. 

Spectral reduction was performed using stan- 
dard irafU tasks. The blue data frames were 
bias subtracted, registered, combined, cosmic ray 
cleaned, distortion corrected and wavelength cal- 
ibrated. Relative flux calibration was performed 
using three standard stars; absolute calibration 
was not attempted since conditions were non- 
photometric and obtaining EWs only requires rel- 
ative calibration. 



2 IRAF is distributed by the National Optical Astronomy 
Observatory (NOAO), which is operated by the Association 
of Universities for Research in Astronomy (AURA), Inc. 
under a cooperative agreement with the National Science 
Foundation (NSF). 



Sky and galaxy background subtractions were 
performed seperately for each of the clusters. Two 
slits (not shown in Figure 1) in our slitmask were 
configured to obtain the sky background; averag- 
ing these and boxcar smoothing them in the spa- 
tial direction gave us a single sky background im- 
age. The sky background image was then spatially 
scaled and subtracted from each of the cluster im- 
ages. Next, the one-dimentional spectra of the 
clusters were extracted with an aperture of about 
1". Two adjacent apertures of the same width 
were simultaneously extracted and combined to 
yield the background galaxy spectrum. Both the 
cluster spectrum and the background galaxy spec- 
trum were smoothed by 9 pixels (— 9.8 A). The 
continuum and the emission spectra of the clusters 
were found to be clearly defined, and the galaxy 
background had no steep gradient for each of the 
clusters, thus allowing a simple linear subtraction. 
Visual inspection suggests that this subtraction ef- 
fectively removed the galaxy background from the 
cluster spectra. Figure 1 shows the spectra of the 
clusters, plotted as relative flux versus observed 
wavelengt h. Note that c luster candidate S3, or 
No. 32 in Ide Grijs et al. ( 20031 ) . is a foreground 
star. 

3. SPECTROSCOPIC RESULTS 

3.1. Ages of Star Clusters in NGC 4676 

The spectra of most clusters show a nearly flat 
continuum with a combination of strong nebu- 
lar emission and underlying stellar absorption fea- 
tures, suggesting populations of very young mas- 
sive stars surrounded by ionized gas. Two of 
the spectra, cluster SI and S9a, show Balmer ab- 
sorption lines along with relatively weak Ca II K 
lines, indicating t hat these clusters are dom inated 
by A- type stars (|Schweizer fc Seitzerlll998h . The 
Balmer jump and the prominent optical contin- 
uum of these two clusters are signs of faded young 
massive stars and older ages. 

The Balmer lines (H/3, H 7 , U5, H8 and H9) 
were measured using the SPLOT and NGAUSS- 
FIT tasks in IRAF. In order to separate the emis- 
sion and the absorption features of some Balmer 
lines, both features were simultaneously fit with 
Gaussian profiles. The equivalent widths (EWs) 
were then calculated and redshift corrected, i.e. 
EW rest = EW b se ™ e d/(l+z). Age estimates for 
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the clusters were derived by comparing the mea 
sure d EW of each line with the synthesized E W 



from l Gonzalez Delgado. Leitherer fc Heckman (1999, 
hearafter GLH99), assuming a Salpetcr IMF and 
solar metallicity. Note that the nebular emission 
of B aimer series were used to derive ages less than 
20 Myr (c.f. Table 7 in GLH99), and for clus- 
ter SI and S9a the absorption features were used 
(c.f. Table 1 in GLH99). Averaging estimates for 
different lines yielded the age and the statistical 
error for each cluster shown in Table 1. 

Although the signal to noise ratio was not high 
enough to detect most He I lines, He IA4471 was 
observable in five of the spectra. The spectra of 
cluster S2, S5 and S10 showed emission features, 
whereas cluster S7 and S8 showed only absorp- 
tion features. The EWs of the He IA4471 line 
were measured using SPLOT, redshift corrected 
and compared with the synesized EWs using the 
same method described above (c.f. Table 3 and Ta- 
ble 7 in GLH99). Age estimates based on this line 
and the calculated uncertainties from propogation 
of errors are also shown in Table 1 . The age esti- 
mates from Balmer lines and He IA4471 line are 
in a good agreement with each other. 

3.2. Metallicity Results 

We adopt the chemical analysis method of 
Kobulnickv fc Kewlevl (|2004l hereafter KK04) to 
determine metallicity. The EW ratio of the colli- 
sionally excited [O Il]A3727, [O III]AA4959, 5007 
emission lines relative to the H/3 recombination 
line, known as R23, allows an estimate of gas- 
phase oxygen abundances, 12 + log(0/H), of the 
young clusters. Instead of traditional flux ratio, 
this method using EW ratios has the advantage of 
being reddening independent to first order. The 
EWs of the [O II] and [O III] lines were measured 
using the SPLOT task in IRAF and were redshift 
corrected. The EWs of H/3 used for the following 
metallicity calculation were the same as those used 
for the age estimates, where the emission features 
were corrected from the absorption features. 

The oxygen abundance, however, is known to 
be degenerate with the R23 diagnostic line ratio, 
giving both a high "upper branch" and a smaller 
"lower branch" metallicity estimate (e.g. Figure 7 
in KK04; the separation between the two branches 
occurs at 12 + log[0/H]~ 8.4). In this work 
we calculated the oxygen abundances for both 



branches using the analytical expressions of Mc- 
Gaugh(1991) in KKO^l To break the degener- 
acy, we used the EW of the [O II] A3727 line as 
the discriminator. L. J. Kewley & H. A. Kobul- 
nicky (2007, in preparation) suggest that the oxy- 
gen abundance would have values represented by 
the upper branch if the EW([0 II]) < 50 A or by 
the lower branch if the EW([0 II]) > 200 A. Note 
that for cluster S10, the EW([0 II]) lies in be- 
tween these limits so that its oxygen abundance 
is still uncertain. In Table 1 we list the oxygen 
abundances of each cluster under the "lower" and 
"upper" columns, followed by the estimated uncer- 
tainties due to the H/3 absorption features. These 
uncertainty estimates do not include the error in- 
troduced by the model uncertainties, typically ~ 
0.1 dex, in the theoretical calibrations. Our best 
guess for the metallicity, depending upon which 
branch it lies on, is listed in the "adopted" col- 
umn of Table 1. 

For most of the clusters, their adopted oxygen 
abundances are greater than 8.91. This value cor- 
responds to a metallicity of 1.6 ZgE]. Cluster S2 
has an adopted oxygen abundance of 7.34, which 
corresponds to a metallicity of 0.04 Zq. These 
are distinctly non-solar metallicities, however our 
age estimates were based upon a model assum- 
ing a solar value. Therefore corrections to the 
derived ages needed to be considered. The ef- 
fect that these metallicities have on our age es- 
timat es is manifested in th e evolution of the H/3 
EW (Leit herer et al. Il999h . For ages less than 



500 Myr, the main-sequence stars that dominate 
the luminosity have lower effective temperature at 
high metallicity; consequently the nebula produces 
weaker H/3 emission lin e s. Co mparing to the re- 
sults of iLeitherer et al" (1999) our age estimates 
for clusters with higher metallicity were oversti- 
mated by approximately 1 Myr. In contrast, the 
low metallicity for cluster S2 produces an under- 
estimate for the age by about 1 Myr. 



3 Equat ion 15; also see iKobulnicky. Kcnnicut t fc Pizagnol 
(1999) Equation 8 for expressions of oxygen abundances 
on the lower branch. 

4 The metal mass fraction, Z, is related to the oxygen a bun- 
dance by Z ~ 29 X ioll2+log(0/H)]-l2^ for the stanc iard 
solar ab undance distribution with t he solar oxygen abun- 
dance of Allcndc Pricto et al. (2001), which yields a solar 
metallicity of Z Q ~ 0.015 and 12 + log(O/H) = 8.72. 
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4. DISCUSSION 



5. CONCLUSION 



We obtain reliable age determinations from our 
spectra of the star clusters in the Mice, and as a 
collateral re sult calculate their metalli c ities f or the 
first time. Kewlev. Geller fc Barton ( 20061 ) sug- 
gests that low nuclear metallicities in close pairs, 
atypical of the metallicity gradients of late-type 
spiral galaxies, may be evidence for tidally in- 
duced gas flows. The galaxy interactions cause 
gas flows toward the central regions, delivering 
less enriched gas from the outskirts of the galaxy. 
Our metallicities show a nearly flat distribution 
along the northern tail of t he Mic e, consistent with 
Kewlev. Geller fc Bartonl (|2006f) scenario. This 
flat metallicity distribution may suggest gas mix- 
ing within the tidal tail. 

Barnes! ( 2004 ) has performed numerical simu- 



lations using two star formation rules which yield 
markedly different predictions for star formation 
in the Mice. A shock-induced prescription pre- 
dicts a widespread burst of star formation trig- 
gered by large-scale shocks during the first passage 
175 ± 25 Myr ago and a second burst of star for- 
mation, mostly concentrated within the central re- 
gions of the galaxies, developing more recently. In 
contrast, a simulation using a density-dependent 
prescription predicts no global burst of star for- 
mation associated with the passage, but rather a 
gradual increase within the galactic centers. From 
our results, most young clusters formed fairly re- 
cently (< 10 Myr), however two were found to have 
ages of ~ 170 Myr, which suggests that they most 
likely formed during the first passage of the Mice. 
These two older objects are located in the tidal 
tails of the pair, which is consistent with the spa- 
tial distribution of star formation predicted by 
shock-induced models. 

Young star clusters formed in mergers may 
evolve into globular clusters, illuminating the 
origin of the high specific freq uency of globu- 
lar clusters in elliptica l galaxies (jSchweizer 1986; 
Ashman fc Zapfl Il992h . Cluster SI appears un- 
resolved on the image, and with an estimated 
initial mass of 10 6 M Q (assuming a Salpeter IMF 



Leitherer et al 



with low mass cut-off of 0.15 M Q ; _ 
1999). the cluster is massive and luminous enough 
to qualify as a young globular cluster 



Here we present the first results on the ages 
and metallicities of young clusters in the Mice. 
These ages are obtained by comparing the EWs 
of the Balmcr lines to model values, and the oxy- 
gen abundances are calculated using standard line 
diagnostics. Most of the cluster spectra show 
strong Balmer emission lines, implying estimated 
ages from 2.5 to 7.1 Myr. Two of the spectra lack 
emission features, yielding estimated ages of 163 
and 172 Myr. These ages are comparable to the 
dynamical age of the tidal tails and are consis- 
tent with shock-induced star formation triggered 
during the first passage of the Mice. The calcu- 
lated oxygen abundances for the clusters along the 
northern tidal tail may indicate gas mixing within 
the tail. Further detailed comparison with the pre- 
dicted age distributions of the star clusters from 
simulations will be presented in a future paper. 

We thank Gregory Wirth for expert assis- 
tance at the telescope, and the anonymous ref- 
eree for useful comments that helped improve the 
manuscript. The authors recognize and acknowl- 
edge the very significant cultural role and rever- 
ence that the summit of Mauna Kea has always 
had within the indigenous Hawaiian community, 
and are grateful to have had the opportunity to 
conduct observations from this mountain. L.-H. 

C. gratefully acknowledge the support from the 
STScI through grant HST-GO-09822.07-A, P.I. 
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Table 1 

Property, Estimated Age, and Metallicity of Star Clusters in NGC4676 





Other 


RA (J2000) 


DEC (J2000) 


F475W a 


Age b 


Age c 




log(0/H) + 12' 


i 


No. 


ID a 


(hh mm ss.ss) 


(dd mm ss.ss) 


(HST mag) 


(Myr) 


(Myr) 


Lower 


Upper 


Adopted 


Sl e 


40 


12 46 12.15 


30 43 06.86 


22.51 


172±33 










S2 


35 


12.33 


25.76 


23.26 


2.5±0.54 


2.6±0.36 


7.34±0.03 


8.98±0.02 


7.34±0.03 


S3 f 


32 


09.73 


21.00 


23.55 












S4 


30 


11.53 


36.00 


22.52 


5.5±0.06 




7.47±0.02 


8.95±0.01 


8.95±0.01 


S4a 


31 


11.57 


34.46 


21.80 


6.0±0.20 




7.56±0.03 


8.91±0.01 


8.91±0.01 


S5 


26 


10.18 


48.62 


19.68 


4.1±0.08 


4.8±1.94 


7.33±0.01 


8.96±0.01 


8.96±0.01 


S6 


24 


10.16 


55.37 


20.15 


5.9±0.09 




7.76±0.17 


8.83±0.11 


8.83±0.11 


S7 


8 


10.51 


30 44 34.75 


20.34 


6.6±0.42 


6.6±3.15 


7.53±0.03 


8.92±0.01 


8.92±0.01 


S8 


6 


10.51 


53.30 


20.66 


7.1±0.10 


7.7±2.39 


7.43±0.16 


8.97±0.07 


8.97±0.07 


S9 


5 


10.57 


30 45 16.42 


21.68 


5.9±0.04 




7.58±0.11 


8.92±0.06 


8.92±0.06 


S9a 




10.62 


14.09 




163±35 










S9b° 




10.65 


12.21 




5.6±0.11 




7.62±0.04 


8.96±0.02 


8.96±0.02 


S10 


2 


10.56 


31.62 


22.97 


2.7±0.22 


2.9±0.27 


7.52±0.10 


8.86±0.04 





a Data from lde Griis et al.| 1120030 . 

b Derived by comparing EWs of Balmer lines with GLH99, assuming solar metallicity. 
c Dcrived by comparing EWs of He IA4471 line with GLH99, assuming solar metallicity. 
d Calculatcd from KK04 and Kobulnicky et al. (1999). 
c No Oxygen nebular lines were observed. 
*A foreground star. 
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Fig. 1. — ACS g-band image of the Mice, with blow-ups of star cluster S9, S9a and S9b on the right and S4 
and S4a on the left. Boxes represent the slits used in the spectroscopy. Age of each star cluster is shown with 
its spectrum, plotted as relative flux vs. observed wavelength (A). Markers in the spectra are the Balmer 
series. Note that S3 is a foreground star. 
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